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Abstract: A method is developed for calculating the jet mass distribution at hadron 
colliders using an expansion about the kinematic threshold. In particular, we consider 
the mass distribution of jets of size R produced in association with a hard photon at the 
J> \ Large Hadron Collider. Expanding around the kinematic threshold, where all the energy 

goes into the jet and the photon, provides a clean factorization formula and allows for the 
resummation of logarithms associated with soft and collinear divergences. All of the large 



logarithms of jet mass are resummed at next-to-leading logarithmic level (NLL), and all 
the global logarithms at next-to-next-to- leading logarithmic level (NNLL). A key step in 



OO 
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CN ■ the derivation is the factorization of the soft function into pieces associated with single 

scales and a remainder which contains non-global structure. This step, which is standard 
in traditional resummation, is implemented in effective field theory which is then used to 

^ " resum the large logarithms using the renormalization group in a systematically improvable 

H '. 

' manner. 



1. Introduction 



The properties of jets produced at the Large Hadron Collider (LHC) have been coming 
into increased focus over the last few years, both from theoretical and experimental points 
of view. Jet properties can be used both to understand QCD and to find and discriminate 
among models of new physics. Indeed, many methods involving jet substructure have 
been developed over the last few years that can be used in new physics searches [THT3]. 
These studies are almost all based on Monte-Carlo simulations. As the precision of the 
experimental measurements improves, precision calculations beyond the level of current 
Monte-Carlo generators will become increasingly important. In this paper we provide a 
framework in which the simplest jet property, jet mass, can be computed in a systematically 
improvable way. 

Precision QCD calculations of jet shapes, such as jet mass, are difficult for a number of 
reasons. Calculations at fixed order in perturbation theory are only useful at large values of 
jet mass in the tail of the distribution. Near the peak, large Sudakov double logarithms be- 
come increasingly important. The resummation of these logarithms is critical in achieving 
distributions with even qualitative agreement with data. However, the resummation of jet 
mass in a realistic hadronic collider environment is challenging due to the large number of 
variables: the jet size R, the jet algorithm, the beam remnants, the initial state radiation, 
hadronization, underlying event, etc. 

There has been flurry of work on jet mass resummation over the last few years |14I - 
126] , In particular, it has been demonstrated that large logarithms of jet masses can be 
understood and in some cases resummed using effective field theory techniques. In this 
paper we consider the simplest QCD event shape, the jet invariant mass, and we focus 
on a very simple event topology: a jet produced in association with a hard photon. This 
process provides a natural setting to generalize similar work that has been done for e + e~ 
colliders. Furthermore, such a simple final state allows us to explore issues such as dealing 
with the beam remnant and the jet size, before attempting to extend these results to more 
complicated multi-jet final states. 

Many of the ingredients necessary for jet mass resummation in direct photon produc- 
tion have already been provided in previous work that calculated the direct photon px and 
rapidity distribution using threshold resummation [27] . For the px spectrum inclusive over 
the jet properties, one can treat the jet as everything-but-the-photon. Since the beam rem- 
nants are included in the jet, the factorization formula in the inclusive case is particularly 
simple. In this paper, we want to calculate a cross section which is differential in the jet 
mass. Thus, we must separate the soft radiation into an in-jet region and an out-of-jet 
region, which leads to a modified factorization formula. 

Although the jet mass distribution in this paper is calculated by expanding around the 
machine threshold limit, where the jet and photon have the maximum possible energy, the 
results will be accurate well away from threshold, at phenomenologically relevant values of 
the jet momentum, due to dynamical threshold enhancement [28 -30J. Dynamical thresh- 
old enhancement works because of the extremely rapid fall-off of the parton distribution 
functions (PDFs) as momentum fractions approach unity. This fall off forces radiation in 
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realistic events to have similar properties to radiation in events near threshold. In particu- 
lar, the large logarithms associated with soft and collinear singularities are determined by 
an effective maximum energy, set by a parton threshold, rather than the hadronic thresh- 
old. The power of threshold resummation has been confirmed by direct comparison with 
data in many cases. Here, in the absence of direct photon jet mass data, we confirm it by 
comparing to PYTHIA [31]. 

When dealing with an exclusive observable such as the jet mass, one is forced to 
contend with multiple relevant scales: the hard scale of the partonic collision, the jet 
scale associated with collinear radiation in the jets, the factorization scale associated with 
collinear radiation in the beams, and various soft scales. Through the use of the soft- 
collinear effective field theory (SCET) [321133] . each of these scales is associated with a 
function that accounts for the physics of each scale. The multiple soft scales are particularly 
troublesome to deal with because the factorization formula coming out of effective theory 
does not distinguish them: SCET gives you only one soft function, with multiple scales. 
Previous work on multi-scale soft functions in SCET has observed that, in the absence 
of an exact two loop soft function calculation, improved agreement with QCD can result 
from assuming a factorized soft function |2ip 25j. In this paper, we make further progress 
in understanding this refactorization by taking inspiration from results on resummation in 
perturbative QCD. 

Previous work on jet mass resummation in perturbative QCD, such as in [23 ; 26,34 
has approximated the jet as having only collinear radiation and soft-collinear radiation. 
These two ingredients combine into an infrared finite object that is a type of jet function, 
but different from the jet function used in SCET. This simple approximation leads to 
good agreement with CDF data in [35] and indicates that the dominant large logarithmic 
corrections are contained in the collinear and soft-collinear sectors. To go beyond the 
leading-logarithmic accuracy of [23U26j . we need to also include initial state radiation and 
color coherence effects. These are naturally part of the threshold calculation we provide 
here. We therefore factorize out of the full soft function a part associated with radiation 
which goes entirely into the jet, which has its own natural scale. We give a precise operator 
definition of this regional soft-function that, in particular, contains all the modes that are 
simultaneously soft and collinear to the jet. Thus our results will reduce to previous results 
at leading-logarithmic level, but are systematically improvable. 

Our systematic improvements over previous work include the full next-to-next-to- 
leading logarithmic (NNLL) resummation of the global logarithms, including initial state 
radiation and color-coherence effects. One difficult unsolved problem in jet observables 
which are not completely inclusive is that of non-global logarithms (NGLs) [36H39] . Al- 
though we have not succeeded in resuming NGLs in this work, and thus our distribution 
is only formally valid at the NLL level, our factorized soft function clarifies the role of the 
missing terms. In an related observable, jet thrust |21U25| . the anomalous dimension of the 
analogy of our regional soft function is proportional to the leading non-global logarithm. 
We use the correspondence to estimate the size of NGLs for jet mass by varying the ap- 
propriate coefficient in our expressions. The effect of non-global logarithms is significant 
in the peak region and must be understood better for further systematic improvements. 
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Figure 1: Event topology of direct photon production. 



Because the factorization formula in SCET gives different objects with different as- 
sociated scales, we can resum large logarithms by matching and running between these 
scales. It is an advantage of the effective theory approach that these scales are manifest 
throughout the calculation. Due to the convolution with the non-perturbative PDFs, scale 
choices cannot be made analytically. We therefore choose our scales numerically, following 
the approach of |30j . As emphasized in |40j and [U], it is an advantage of the effective 
theory approach that natural scales automatically appear through this numerical proce- 
dure, in contrast to fixed order calculations where they can only be guessed. Our results 
are compared to the output of pythia [12], with very good agreement. 

2. Kinematics and the Observable 

In this section we briefly review some kinematics of direct photon production in a hadron 
collider and set up our notation. We will follow the conventions in [27j . and we summarize 
the most important kinematic relationships. Once the kinematic variables are defined, we 
will define our jet definition and the specific observable that we calculate in this paper. 

2.1 Hadronic and Partonic Kinematics 

We will denote the momenta of the incoming hadrons as and P^\ an d the momentum of 
the photon as p^. The observable will be categorized by the properties of the hard photon, 
namely its transverse momentum and rapidity y. Formally, the results of this paper 
will be most accurate when p? is near the machine threshold limit, 



where Eqm = \J (Pi + P2) 2 . A majority of events will have photon transverse momenta 
that are far below the machine threshold limit. However, as we will see, the calculation 
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remains valid for transverse momenta much less than machine threshold due a phenomenon 
called dynamical threshold enhancement, discussed in Section 

A schematic showing the event topology for the collision is given in Figure [TJ The hard 
process involves the collision between two initial state partons, which have momentum 
Pi = XiP- 1 , where the Xi is the momentum fraction of the constituent parton and the 
subscript denotes the hadron from which it came. There are two channels that are relevant 
at leading order in a s , the annihilation channel, qq — > 75 and the Compton channel, 
qg — > "fq. It is convenient to use the variables w and v, given by 

1 PT y PT -y , n ^ 

xi = - e y , x 2 = -e y . 2.2) 

w Eqmv Ecm{1-v) 

The variable v is related to the scattering angle between the photon and the beam, 9*, in 
the partonic center of mass frame through 

v = -(l + cos0*). (2.3) 

The other variable w characterizes how much the (partonic) event deviates from exact 
2—7-2 kinematics due to additional QCD radiation in the initial or final state. When 
w = 1, there is no additional radiation and the partonic final state consists of a photon 
and either a quark in the case of qg — > qj or a gluon in the case of qq — > g^j. 

When discussing the threshold limits, we will discuss the (total) hadronic invariant 
mass 

Mi = (P 1 + J P 2 -p 7 ) 2 , (2.4) 
and the (total) partonic invariant mass 

m\ = {x 1 P 1 +x 2 P 2 -p 1 f . (2.5) 

These give the invariant mass of all particles in the final state with the photon's momentum 
removed from the hadronic or partonic final state, respectively. The hadronic variable Mx 
includes the beam remnant, whereas the partonic version does not have any information 
about the beam, and therefore < mx < Mx- In terms of the momentum fractions of the 
partons and the photon transverse momentum and rapidity, the invariant masses are 

M 2 X = E 2 CM - 2p T E CM cosh y = E 2 CM (l - -J^J (2.6) 

and 

9 p% 1 — w , _ s 

2 = _^T / 2 7 x 

X 1-V W V ' 

In the partonic threshold limit, w — >■ 1, there is no phase space volume for the initial or 
final state radiation, and so the the final state consists of a single scattered parton and a 
photon. After removing the photon, only a massless parton remains and so the partonic 
mass vanishes in this limit (mx —> 0). The partonic cross section is singular in this limit, 
and the resulting singularities take the form of In 2 mx ~ ln 2 (l — w). In the next section, 
we will introduce a factorization theorem for the partonic cross section, that will resum 
the logarithms of mx that appear as a™ln m (l — w) in the differential cross section. 
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2.2 The observable 

The observable we consider is the invariant mass of the hardest jet in events with a high 
px photon. Rather than use one of the standard jet algorithms to define the jet mass, we 
use a simple definition in order to make the analytical calculations more tractable. The jet 
mass used in this paper is defined through the following procedure. 

1. Remove the photon from the event record. 

2. Cluster the event using the anti-fc^ jet algorithm. 

3. Select the jet with the hardest px, and find its 3-momentum pj. This defines the jet 

axis as hj = r^r. 

J \vj\ 

4. Define m,R to be the invariant mass of all of the radiation that lies within a cone of 
half- angle R centered on hj. 

The observable we consider is t^riji ' within some range of px and rapidity for the jet and 
photon. Our calculation is not sensitive to the difference in pt between the jet and the 
photon, so we put a cut only on the photon px, integrating inclusively over the px of the 
jet. For concreteness our numerical results will have central photons (y = 0) with a given 
px, integrated over jet rapidities between —1 and 1. 

Although our R differs from a usual collider R, which is defined in terms of rapidities 
and azimuthal angles, the difference is small for central jets that dominate at large px- Our 
method of calculation, and our factorization theorem, would also work for more standard 
jet algorithms, such as the kx algorithm. However, most standard algorithms suffer from 
clustering effects |244l43fl4*6] which are not completely understood. Extending the results 
of this paper to more standard jet algorithms will be explored in future work. 

To calculate the jet mass accurately in QCD, one would like the relevant events to 
come from the hard process we are considering, namely a single photon plus a single 
parton. However, generically, there can be contributions from events with a hard photon 
and 2 or 3 jets. The simplest way to ensure single jet kinematics is to use a px veto on the 
second hardest jet. This is easy to do experimentally or when using an event generator, 
but in QCD introducing a new scale makes the theoretical calculation significantly more 
complicated. In particular, non-global logarithms of the jet mass relative to the veto 
scale will be introduced. Although these non-global logarithms might be a small effect 
numerically, an advantage of the threshold resummation is that one does not have to 
introduce an extra scale. Instead, the scale is implicit in the out-of-jet soft radiation which 
is integrated over. 

By expanding around the threshold limit and demanding large px for the photon, we 
can ensure that we are in the photon-plus-single-jet configuration without introducing a 
new veto scale. When the photon's px is on the order of the hard scattering scale, the 
radiation outside of the hardest jet must be soft, which prevents the formation of a second 
hard jet. By suitably restricting the photon px, the events will contain a single recoiling jet 
and we are safely in the regime where the factorization theorem is valid. To put it another 
way, the fact that dynamical threshold enhancement has been shown to work for the direct 
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Figure 2: The photon's momentum and rapidity is denoted by p 7 and y. respectively. The jet 
cone of half angle R, centered around the direction of the jet nj = (l,hj). The restrictions on 
the jet rapidity yj is shown by the dashed gray line. Soft radiation inside the jet contributes to 
ha = n J ' -fx'"! whereas the radiation outside contributes to k Q ut = nj ■ Px°ut 



photon pt distribution in [27] indicates that the single jet configuration dominates. To be 
clear, the threshold expansion does not obviate the need to deal with non-global logarithms, 
as you would with a jet veto, but it does avoid having to deal with an explicit scale. 

At the partonic threshold limit, where the process only has one parton and a photon in 
the final state, the pt and y of the photon constrains the scattering angle of the outgoing 
jet. To see this, consider momentum conservation along the direction of the beam, taken 
to be the z direction. 

PJz = (xiPi + x 2 P 2 - Pl ) z (2.8) 

After using the definition of rapidity, y = | In ^ p p z , we can relate the rapidity of the jet 
yj to the pt and y of the photon and to x\ 2 via 

sinhyj = ^^-(xi - x 2 ) - sinhy , (2.9) 
2p T 

Since < x\p < 1, yj must lie within a range given by 

CM _ gjjj^ ^ sinhyj < — — sinhy (2.10) 

2p T 2p T 

In order to perform threshold resummation, we must be careful to keep yj small enough 
to ensure the jet does not enclose the beams. This restriction is sinhyj < cot R, and so if 
we take \yj\ < 1, this implies that R < 0.7. 



-6- 




Figure 3: The integration region 1Z is shown by the red lines. When the jet rapidity is restricted, 
the region is further constrained to be inside the blue lines, shown here for y cut = 1. The dashed 
green line outlines the region that most strongly contributes. This is demonstrated by the right 
hand plot, which shows the rapid fall off of x\f u /p{xi)X2fu/p{ x 2)- 



3. Differential Cross Sections and Factorization Theorem 

The differential cross section for the process N\ + N 2 — > 7 + X can be written in the form 

d PT dydml = ^ JJ n dvdw Wa/* A*)] M/N a M\ > (3-1) 

where iVj are the incoming partons, f a /Ni( x ii A*) is the- probability distribution function for 
finding parton a in hadron iV, with momentum fraction Xi, and the sum is over the different 
partonic channels. The region of integration 7Z in the w — v plane is given by 

1 + ™MA < v < 1 _ _PT_ e -v 
f^e-y + m 2 R /p T E CM 



PT 



°y Pt ^ P 2 T 

< w < 



v E C m " Pt + (1 - v)m 2 R 



e 



!J 



— — — < id < — 1 — — . (3.2) 

u(2smhy + 2smhy cut + j-^e f) u(2smhy - 2smhy cut + — e v) 

The first two constraints account for < Xi < 1 and m|j < , which is the region 
enclosed by the red lines in Figure [3j the third constraint comes from a finite rapidity 
interval for the jet: — y cut < y,j < y C ut, drawn as the blue lines in Figure [3l The expression 
d2 T— 3- gives the partonic cross section for a + b — > 7 + X. For a fixed w and u, the 



duid^dm^j 



partonic cross section, at NLO, will diverge logarithmically as rriR — > 0, and in this regime 
these logarithms can invalidate an expansion in a s . 

The reason that the threshold expansion works away from threshold is that the parton 
distribution functions fall rapidly as the momentum fractions approach unity (xj —> 1). In 
Eq. ()3.5p . the partonic differential cross section is weighted by the product of the PDFs. 
As a representative example, the right side of Figure [3] shows a plot of the product 

xif u /p(xi, Hf)x 2 f u /p{x2^f) , (3.3) 
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which would be necessary for process uu — > 75. In the plot, we have chosen Ecm = 8 TeV, 
Pt = 2 TeV, y = 0, and [if = \/H = ■ v /xiX2-Ecm- The plot clearly shows the extremely 
rapid fall off of the PDFs, and this demonstrates that, for a fixed value of rriR, the region 
in w — v space is weighted far more heavily along the line given by mx = mR. 
The hadronic invariant mass can be written as 

2 

M\ = ^ + E 2 CM [(1 - x x )v + (1 - x 2 )(l - «)] . (3.4) 

When Mx —> 0, this forces both mx —> and x\ 2 - ^ 1; however, the fall off the PDFs 
favor small mx even when Xi is not close to 1. Since mR < mx, this also means that 
small mR is favored. The partonic cross section is singular in mjj, which manifests itself as 
powers of In mji, and for small mR, these logarithms are large and must be resummed. This 
enhancement of the singular region for the partonic cross section means that resummation 
of mR is important, even when we are not near the machine threshold. This effect is known 
as dynamic threshold enhancement [28H30] . 

3.1 Factorization of the partonic cross section 

When the jet invariant mass is small compared to the pr, the jet is highly collimated and 
the process can be described within the framework of the Soft-Collinear Effective Theory 
(SCET) [32,33,4TJI38]. Rather than give a formal derivation of the factorization formula 
as in |27j . we will focus on how the factorization formula must be modified to become 
differential in the jet mass. 

The partonic differential cross section can be calculated via 

H^/TT f 2E c\Jc 771 ^ 777 ^ 

dwdvdm 2 R = *> H(PT,v,») J o J — J{m 2 R -{2Ej)k,») S(k, ^ ■ (3-5) 

The three functions H, J and S are the hard, jet and soft functions, which account for the 
physics associated with the different scales involved. The quantity a(v) is defined through 
the leading order QCD result via, 

dwdvdm 2 = dab ^ S ^ m ^ ~ m ^ ' ^ 3 ' 6 ^ 

where the tree level parton cross sections are given by 

<7qq[ v ) = nu eia e q a s (ii)-^- I — 

1 /. n9n V 



o qg {v) = ira cm e 2 q a s (n) — (1 + (1 - v) 2 )^-^ . (3.7) 



The hard function is given by the square of the Wilson coefficients that arise when matching 
the QCD amplitude onto the effective field theory operators. It describes the short distance 
physics at an energy scale fi^ ~ pr- The jet function comes from integrating energetic 
modes that are collinear to the jet direction that have virtuality comparable to mR. In the 
threshold limit mR <C PtR, so the collinear radiation is insensitive to the jet boundary. 
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Thus we can use the same inclusive jet function as that in [27]. Finally, once the collinear 
modes are integrated out, we are left with the soft modes, which are described by soft 
Wilson lines. The soft function relevant for the jet mass calculation is different from the 
one used in [27). As we will see, it depends separately on the radiation inside and outside 
of the jet, and also on the jet algorithm. 

The form of the factorization theorem can be understood from power counting in the 
threshold region. In the partonic threshold limit, corresponding to the limit mx — > 0, 
radiation is constrained either to be collinear to recoiling parton, which we call pj, or to 
be soft, which we write as k^. In terms of these momentum regions, the partonic invariant 
mass takes the form 

m\ = {pj + kf = m 2 c + 2k-pj + --- (3.8) 

up to terms of order k 2 <C m 2 where m 2 = p 2 is the collinear mass, whose distribution is 
given by the jet function. 

The jet invariant mass m R only contains the part of radiation inside the cone. Since 
all the collinear radiation is in the cone, we have only to split the soft radiation, which we 
write as k^ = k^ n + fc„ ut , with k^ n in the cone and k^ ut out of the cone. Then, 

m R = (PJ + kin) 2 = m 2 c + 2fcin ■ Pj ■ (3.9) 

and so 

m x = m R + 2 ^out • Pj • (3.10) 
Now, the jet momentum pj accounts for radiation that is collinear to recoiling parton, and 
so it can be written as 

Pj = Ejiij + residual . (3-H) 

where rij = (1, hj) is a lightlike vector that is directed along the jet direction. The residual 
momentum gives a power suppressed contribution to the jet mass compared to the leading 
component, and thus can be dropped. The expressions for mx and m R then become 

m x = m R + 2Ej(nj ■ fe out ) 

m 2 R = m 2 c + 2Ej(nj ■ k- m ) . (3.12) 

Therefore, the only component that enters the soft function is k = nj-k. These expressions 
explain the form of the factorization theorem, Eq. (|3.5p . which can be written as 

d 2 & f 
——— — r = w a(v) H(p T ,v,n) / dfc in d/c out dm 2 J(m 2 )5'(/ci n ,fc ou t) 
awavam R J 

x5(m 2 R - m 2 - 2Ejk in )5{m 2 x - m| - 2Ejk out ) . (3.13) 

Furthermore, we see why we need a soft function that depends on the projection of the soft 
momentum on the jet direction separated into in-cone and out-of cone components, which 
we denote k{ n and k OU f 

The precise definition of the soft function we need is, for annihilation channel, 



s qq (k in ,k out ) = ; ^^^|(x s |T[y 1 V J t%V 2 (o)]^^ , 

(3.14) 



C F N C 
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and for Compton channel 
1 



Sqg ( fcin j &out ) 



C F N C 



x s 



^2\{X S \T Y?Y 2 t a Y}Yj(0) \0)\ 2 5(nj-P x ^-k in )5(nj-P Xs » t -k out ) 



(3.15) 

where Pjf in an d Px° ut denote the total momentum of the soft radiation propagating inside 
and outside the jet, respectively. Also, Yi is a soft Wilson line directed along the 7ij = (1, hi) 
direction: 



Y^{x) = Pexp fig J ds m ■ A s (riiS + 



x) 



(3.16) 



where i = 1,2, J and P denotes path ordering and A s = A^T a are soft gauge fields in the 
fundamental representation. 

3.2 One-loop soft function 

We now calculate the required soft function to order a s . The diagrams needed for the 
calculation are represented in Figure HI The virtual diagrams are not shown as they give a 
scaleless contribution when using dimensional regularization to regulate both the ultraviolet 
and infra-red divergences. The non- vanishing diagrams for the qq — > gj channel give 



Sqq(kjjx, &out) M) 



Cf - -xCa 



8ira s (ix e 



(27T) 



d-1 



47T 



d d q 5(q 2 )9(q )M(k in ,k OVLt ;q) 



ni ■ n 2 + 1 nj ■ ni + 1 n j • re 2 

(ni • g)(n 2 • <?) 2 {nj-q)(ni-q) 2 (nj ■ q){n 2 ■ q) 



,(3.17) 



where q + = nj ■ q, and q = nj ■ q with nj = (1, —hj), and the measurement function is 
given by 



M(k in ,k ont ;q) = & 



q~ + <T~ 
+9 ( cos P 



cosP)5(fc in - g + )<5(/cout) 

£(fcin)<5(fc ou t - 



(3.18) 
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Similarly, for the qg —> q"f channel we have 



Sqg(k- m , k ou t, /i) 



8ira q 



<lc A - 



(2Tr) d ~ 1 V 4vr 
ni ■ n 2 /„ ! 



A % 2 )%o)-M(MnAut;g) 



+ 



(c,--c x ). 



+ \c A - 



nj ■ n 2 



.(3.19) 

L2~^(ni • g)(n 2 • g) ' 2~*J(n J -q)(n 1 -q) ' 2~ A (nj-q)(n 2 -q)-' 

The soft function will depend on the angle i? and the rapidity of the jet yj. For convenience 
we define 

R 



tan ■ 



P = exp(-yj) 



(3.20) 



With these definitions, the soft functions associated with the two channels at order a s are 



S q q{k- m ,k out ,n) = 6(k in )5(k out ) + 2(^-^6 (k hl ) 5 (k out )(c F - -Ca) In' 



+ 2 
+ 2 



S(k, 



out / 



k'n 

k 



out 



2C A ln 
4C F In ( 



fein l+/3 S 

kout 1 + /3 2 



/i 



/3 



- c, 



<?<? 



}■ 



(1 + /3 2 ) 2 vr 2 
/3 2 6 



(3.21) 



and 



S qg {k hl ,k out ,ii) = S(k in )S(k out ) + 2[^j{5(k iQ )5(k ou t)(-C A ^ In 2 



(1 + /3 2 ) 2 vr 2 - 
/3 2 6 



+ 2 
+ 2 



5(A: 



2C F ln 



where 



Cqq = -C F ln 



-qq 



"an 
<?(fcin) 
&out 

(/3 2 -r 2 )(l -r 2 /3 2 ) 

2^4 



fcin l+/3 S 



</// 



(2C F + 2^)ln(^l±^)-c w ]} 



/3 2 r 



/1 



-CUln 



/3 
/? 2 



(3.22) 



/3 2 - r 2 1 - r 2 /3 2 

Cf In — 5— Ca In 



(1 + /3 2 ) 2 
C F ln 



+ (2C F - C A ) In 



1 + r 2 



/3 2 1 + r 2 

+ C A ln^^. (3.23) 



(1 + /3 2 ) 2 



r 2^2 r 2^2 

The logarithmic divergences involving r and /3 arise when the jet cone encloses the collinear 
radiations from one of the incoming beams. 

Since the factorization theorem contains convolutions between the jet and soft func- 
tions, it is easier to perform the calculation in Laplace space. The Laplace transformed 
soft function is defined as 

k'm \ ( k ou i 



5(Ki n , K ut) A*) 



/ dk hl / dk out exp ( ^— ) exp ( 

Jo Jo v K in eTE/ V 



^jS{ki n , k ou t, fx) ■ 

(3.24) 



Our result is then 

~ S qq(L u L 2 , /i) = 1 + (g) [4(2C F L 2 - CaL\) + 4(L X - L 2 )c w - + (<7 F - ^ A )vr 2 
S w (Li,L 2 ,/i) = 1+ (g) [4((C F + C A )L 2 -C F L 2 ) + 4(L 1 -L 2 )c w + (^Ca)^ 



(3.25) 



- 11 - 



where L\ = In {^j-^j and L 2 = In 



1+/3 2 /tout 

The renormalization group equations satisfied by the soft function in the two channels 



arc 



ds 



qq 



din /i 



In 



1 + /3 2 «i, 



ds 
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dln/j, 



cusp I n 

P V 

1 + /3 2 K in 



4Cf7cus P In 



1 + /? 2 Kout 



2^ Sqq 



2Cf7 cusp l n 



A* 



2(C F + C A )j cnsp In 



/3 // 

1 + /3 2 K out 



'99 • 

(3.26) 



sySqq ^ 



T 



are given to order a 2 in |27] and are known to 



where the expressions for 7 CU sp> 
order ctf. 

To check that the factorization theorem in Eq. (|3.5|) is correct, we have confirmed 
that the //-dependence of the hard, jet, soft functions and PDFs near x = 1 cancel when 
combined. We also compare the singular part of the full QCD distribution at high px which 
we compute using the numerical package mcfm to the distribution produced by combining 
the hard, jet and soft functions at fixed order. This is shown in Figure El 

Although the factorization theorem is correct, it does not automatically guarantee that 
we can resum all of the logarithms of uir. Recall that the soft function in the factorization 
theorem in [27] was fully inclusive and had dependence only on the projection of the soft 
momentum in the jet direction. For the inclusive photon p F spectrum, such a one-scale soft 
function is natural, since the observable only depends on a single scale, pr- In contrast, 
the soft functions given in Eqs. (|3.14|) and (|3.15p depend on two physical scales k[ n ,k ou t 
as well as on the renormalization scale fi. With multiple soft scales, there may not be a 
natural choice for fj,, and therefore the renormalization group evolution may not resum all 
of the large logarithms. Indeed, in a analogous jet mass calculation for an e + e~ collider, it 
was demonstrated first numerically in |21j and then analytically in [25] that the singularity 
structure of QCD was not completely reproduced by expanding the resummed result for 
jet thrust. One expects the same difficulty here. 

The problem is that our soft function anomalous dimension does not depend on jet 
size R. This is simply because the hard and jet functions as well as the PDFs are the 
same as in the inclusive case and do not have a jet size dependence. However, there are 
large logarithms in full QCD that could only be predicted from an i?-dependent anomalous 
dimension. To make this point more concrete, consider the singular part of the cross section 
calculated to 0{a s ) using SCET in the annihilation channel: 



do 



qq 



dwdvdm 2 R 



wo, 



qq 



[v)5(m 2 x 



m R ) 



-2C A T In + c w -(r, 0)T O + T J Q 3 In m 



4tt 



m 



R 



(3.27) 

Since SCET contains all of the singular contributions from full QCD, this will agree with 
the full one-loop QCD result expanded at small mR. However, we should also be able to 
predict these terms only using the anomalous dimensions, since they should be part of 

Ft 

the resummation at NLL level. Since the coefficient of — j- depends on r = tan R/2 through 

m R 

the Cgq(r, (5) term, this is impossible for i?-independent anomalous dimensions. Thus, 
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without further insight, resummation in SCET will undoubtedly miss these logarithmic 
corrections. 



4. Refactorization of the Soft Function 



The soft function depends on the two scales k m and k out . By explicit calculation we found 
it has logarithmic dependence on both of the scales divided by \i already at one loop. 
Since there is not a single scale associated with the soft function, we should not expect 
that a single choice of fi will give the soft function a controlled perturbative expansion. In 
addition, at two loops, one generically expects a complicated dependence on ki n /k out in 
the soft function, as was found in [25j[39], representative of the non-global structure. The 
problem, from an effective theory point of view, is that there are two modes with different 
characteristic scaling described by the same object. Ideally, we could factorize the soft 
function into separate objects which depend only on k- m / /x and fe ou t//x. However, it does 
not appear that the soft function exactly factorizes in that way. 

In the absence of a full refactorization, we will proceed in a conservative fashion fol- 
lowing methods of traditional resummation }50H55| . To that end, we can at least isolate 
part of the soft function whose scaling we already know, namely the part associated with 
radiation that goes entirely into the jet, which will at least contain all of the soft-collinear 
divergences. To do so, we will define an auxiliary soft function that only depends on ki n , 
reminiscent of eikonal jet functions in [50]. There are many ways to do this. A simple choice 
is to define the auxiliary soft function the same way we define the 2-scale soft function, but 
inclusive over the out-of-jet radiation. That is, we take 



£*qq\kin) 



1 



C F N C 



\(X s \T\Y?Yjt a Y}Y 2 (0)] \0}\ 2 5(k in - nj ■ P x , 



(4.1) 



for the annihilation channel, and similarly for the Compton channel. We call this a re- 
gional soft function. The regional soft function was deliberately chosen to have the same 
color structure as the original soft function. 

At one-loop the regional soft function is given by Eqs. (|3.17p and (|3. 19|) with the 
measurement function replaced by 



M(k in ;q) = e( Q q , - cos R)5(k in -q + ) . 
\q + q + 



At order a«, we find 



Sqq(k in ,ri = 5{h n )+2(^[5{k ln )^- l -C A 



hr 



1 

V 

n in 



",<l 



Sqgik^fl) = 5(k in )+2^)\ 



1 



C F 



In" 



/3 2 6 

(1+/3 2 ) 2 7T 2 

P 2 6 



+ 2- 



1 



2C F hx[-—) 



mi 



(4.2) 



(4.3) 
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where c qq - and c qg are given in Eq. (|3.23[) . The anomalous dimensions at order a s , in 
Laplace space, are 



d _ , . _ as 
dln P i Sqq{Kin,tJ ' ) ~ 4vr 



r/ 



din fi 



47T - 



2C A r ln 
2C F r In 



1 + /3 2 ^in 

/3 /i 



4c 3 g(r,/3) s q q(K in ,n) 

4cqg(r,/3) S qg (K hl ,fl) 



(4.4) 



The dependence of the soft function on /c out is still encoded in the full soft function, 
so we define the residual soft function to be what remains after dividing the full soft 
function by the regional soft function. 



Sr(kim ^out? A 4 ) 



S(kn,(J>) 



(4.5) 



At order a s , the contribution from radiation going inside the jet is identical between the 
full soft function and the regional soft function. When the radiation goes outside of the jet, 
it contributes to fc ou t in the full contribution, however its contribution to the regional soft 
function vanishes since the integral is scaleless. Thus the residual soft function at order a s 
only has k out dependence: 

(1 + /3 2 ) 2 - 2 - 



Srqq(hn, kout, fJ>) = S(k ont ) + 2 (^-j |^(fc ut) (pF 



In" 



ft 2 



7T 

If 



4C F ln 



>-out 



ft 



} 

Srqg(k- m ,k ou t,fi) = S (k ont ) + 2 (J^j ^5 (k ont ) (^C F + A 



+ +2 



kn 



(2C F + 2C A ) In 



fcout 1 + ft' 

ft 



In' 



'HI 



(1 + /3 



2\2 



ft 2 



7T~ 

~6~ 



(4.6) 



Thus, the one-loop anomalous dimensions of the residual soft function depends only on the 
Laplace conjugate variable for k out , K Q ut : 



d 



dlnfi 
d 

din// 



'rqij( K out ) A 1 ) — ^ 



4C F r In 



1 + ft 2 K out 



ft 



'rqg 



(«out, A 4 ) 



4tt . 



2(C F + C A )T ln 



1 + ft 2 ^out 



+ Ac qq (r, ft) S r qq(n out , fi) 



ft 



+ 4c qg (r,ft) S rqg (K OVLt ,fi) , 

(4.7) 



Combining these results, we see that we have actually factorized the soft function at 
one-loop, since the regional and residual soft functions each depend on only one scale. At 
two- loops and higher, the regional soft function will still only depend on k[ n , however the 
residual soft function will have both k{ n and A; ou t dependence. Thus, at least there should 
be a natural scale which minimizes the large logarithms in S(k[ Q ,/j>). In particular, this 
scale should correspond to the "see-saw" scale pL s = lAj [ih for the soft radiation associated 



with a particular collinear direction [56] . Such a scale assignment is implicit in much of 
the work on traditional resummation, and here we have made it explicit using SCET. 
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For this refactorization to make any difference, we have to allow the residual soft 
function scale to be different from the scale for the regional soft function. So we write 

S(h n ,k ollt ,fi) = fan) S{k in , fi in )H r (iJ,, fj, out )S r (k- m , feout,Mout) ( 4 -8) 

where the evolution kernels tl(fi,/j,i n ) and II r (/x, // out ) evolve their respective soft functions 
from (jl to separate scales. In Section [6J we investigate how choosing two soft scales affects 
the jet mass distribution by comparing the calculation with and without such a scale 
separation (see Figure [9]) . 

Unfortunately, since the residual soft function will depend on both k\ a and k out beyond 
one-loop, our refactorization is only approximate. At two-loops, the h m /k on t dependence 
of the residual soft function is due to non-global structure. When the scales are widely 
separated, k m <C fc ou t or vice verse, the non-global structure takes the form of a logarithm of 
the ratio of these scales. There have been many studies of non-global logarithms (NGLs) 
[T91IM 1251 1361 ^1291 157] . We will not attempt to calculate the non-global structure for 
this observable, which would require a two-loop calculation, instead we will estimate their 
effects and discuss their phenomenological significance in Section [71 

4.1 Comparing with pQCD result 

Calculating the jet mass distribution at hadron collider has been done by using traditional 
pQCD resummation technique [23j[26l[58]. It's interesting to compare the pQCD results 
and the SCET results obtained in this paper. In Refs. [23, 26], the jet mass distribution is 
described by a process independent quark or gluon jet function, depending on the partonic 
origin of the jet. These pQCD jet functions are calculated to 0{a s ) explicitly, and resumed 
to all orders at NLL level by solving an evolution equation in Mellin space. In their 
approach, both the hard-collinear mode and soft mode contribute to the jet function, in 
contrast to the SCET approach, where jet function is defined with zero-bin subtraction [59] • 
Furthermore, the jet function defined in Refs. [23026] is independent of underlying partonic 
process, and the corresponding jet mass distribution misses the contribution from soft large 
angle radiation between different colored partons. Therefore a simple connection between 
their jet function and our SCET results can not be made. 

To make progress comparing the different formalisms, it is helpful to consider the 
R — >■ limit. In this limit, a process independent combination of jet and soft function can 
be defined in SCET. For example, for the annihilation channel, we have 

Jg CET ( m R> AO = lim / dp 2 dk m J g {p 2 ,v)S g q{k m , fi)5(m 2 B - p 2 - 2Ejk in ), (4.9) 
a R— >0 J 

and similarly for the Compton channel. The fact that in the R — > limit S q g depends only 
on the parton that initiates the jet can be seen explicitly by taking the limit in Eq. ([4. 3D . 
In this limit, only the terms proportional to lnr survive, each of which is proportional the 
jet's Casimir - Cp for quark jet and Ca for a gluon jet. Alternatively, it can be understood 
by the rescaling argument from [25] . To be specific, by rescaling R to 7r/2, the cone of 
the nj-jet becomes a hemisphere. At the same time, all other Wilson lines, except the one 
along the jet direction, are forced to point in the direction of — hj. Therefore, a soft dipole 
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contribution from partons i and J, where % denotes an initial state parton, is identical to 
the hemisphere soft function with the color factor replaced by — 2$ • Tj in the color 
generator notation |6Uj . while dipoles contributes from two initial state partons i and j 
vanish. Since (— T\ ■ Tj — T2 • Tj) = Tj, summation over different soft dipole gives a result 
which depends only on the jet direction and its color content. This explains why the jet 
function constructed in Eq. (|4.9p depends only on the parton that initiates the jet. Indeed, 
at 0(a s ), for quark and gluon jet, we have 

jsc ET{ml „ _ ^ (c A log _ n Cj + i„ /T \ (4 10) 

respectively, where raj is the number of light flavor and Tp = 1/2 in QCD. Eq. (|4.10p 
agrees with the corresponding 0(a s ) jet function in Eqs. (A3) and (A4) of Refs. [26] . 
after power expanding in jet mass, except the term proportional to njTp, which is missing 
there. From the discussion above, we see that the results presented in Refs. |23tl26| are 
not complete at NLL level, missing the contribution from large angle soft cross talk, as 
well as the term proportional to rifTp in the gluon jet function. These extra global NLL 
contributions are taken into account in a more recent work [58], which also includes a 
numerical exponentiation of part of the non-global structure. 

5. Scale Choices 

While the resummed result is formally independent of the scales /i^, fij, /ii n and fi ou t as 
well as the factorization scale we call fi, there is residual higher-order dependence on these 
scales if the perturbative expansions of the hard, jet and soft functions are truncated at 
a finite order. To get a well behaved expansion with minimum hard, jet and soft scale 
dependence, we want to evaluate each contribution at its natural scale where each function 
does not involve large perturbative logarithms. We evolve the hard, jet and soft functions 
from their natural scales to the factorization scale \x as shown in Figure [SJ 
In the resummed distribution, the following ratios of scales appear 

PT m R m X- m R ( r U 

2 , 2 , , . K - 1 ) 

Unlike the fully inclusive case studied in [27], it would be impossible to minimize logarithms 



of these ratios with only 3 scales /i^, /ij and fi s = /Xj n = /i ou t- Instead, we will allow [i- m 
and pi out to be independent, as explained in Section [H When we choose a single soft scale, 
the resummed SCET distribution is hopelessly different from the pythia output - this is 
shown in Figure [9j 

The naive scale choices are 

2 22 

Vh=PT , = m R , /j, in = , n out = ^ , (5.2) 

PT Pt 



- 16 - 



H J S mt S,, 

[i-h , 



Mi 1 

Ms,,!,; 

f^Sin 1 

fi 1 i 1 — y— 

Figure 5: The renormalization group evolution of the hard, jet and soft functions, as shown 
schematically. Each function is calculated at fixed order at its natural scale, and they are evolved 
to a common scale fi. Two scales emerge in the soft sector which results in the refactorization of 
the soft function. 
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Figure 6: The NLO corrections of the hard, jet and soft functions to the jet mass distribution 
as a function of the renormalization scale for mR = 50 GeV and pr — 500 GeV . The natural scale 
of each mode is indicated by the of the extremum of the appropriate curve. The hard and soft out 
scales do not depend on the mass of the jet; they are determined by varying the scales with the 
NLO delta function contributions. The jet and soft; n scales are determined by varying the scales 
with the NLO jet mass distribution. The jet and softi n curves are scaled up to plot with the ones 
for the hard and soft out . 



These would eliminate the large logarithms at each energy scale at the partonic level. 
However, when the partonic cross section is integrated against the PDFs, /U ou t will get 
arbitrarily close to zero and hit the Landau pole singularity of the running coupling. The 
problem is that the natural /x ou t scale choice depends on an unphysical quantity, namely 
m 2 x — m 2 R . We can avoid the Landau pole and minimize the effects of the large logarithms 
by instead choosing the scales numerically to depend only on observables [30]. In fact, it 
is not necessary to eliminate the logarithms in the unphysical, partonic cross section, but 
rather we want to eliminate the large logarithm in the physical, hadronic cross section. 
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20 40 60 80 100 



m R (GeV) 

Figure 7: The jet and softj n scales as a function of ran for px = 500 GeV with two different cone 
sizes. Dots are extrema and curves are our fits. Note that /j,- m is cone size dependent and increases 
when the cone gets smaller. However, /i|//x; n is ran independent, which is a manifestation of the 
factorization theorem. The numerically determined jet and softj n scales fit nicely with power law 
curves (m R in GeV): \i 3 = 1.6 ra R A7 , fj, in (R = 0.3) = ^/2300, fi in (R = 0.5) = /^/6700, and they 
are considerably higher than the naive choices of Hj = ran and /x; n = ra 2 R /pT- 



To determine the natural RG scales, we follow the approach used in [27] [30 |H0| [¥T] . 
We include separately the NLO corrections for the hard, jet or soft functions and vary the 
relevant scale to find out which scale minimizes the variation of the NLO distribution. For 
example, Figure [6] shows these NLO corrections for the 8 TeV LHC with apr = 500 GeV 
photon and m,R = 50 GeV. The extrema of these curves indicate a natural scale for each 
mode, and one can see that there is a natural hierarchy of the various matching scales. 
The jet and softi n scales are jet mass dependent, whereas the hard and soft out scales are 
jet mass independent (their NLO contributions have no dependence on thr). We fit the 
jet mass dependence by power law curves, examples of which are shown in Figure [7] for 
different jet sizes. For example, with R = 0.5, the scale choice we extract are 

fj, h = p T ; , fij = 1.6 m R A7 , fj,- m = , //out = 280 GeV. (5.3) 

We choose the factorization scale to match the hard scale, \l = (a high factorization 
scale is not natural in SCET, however it corresponds most closely to the scales where the 
PDFs have been fit). The renormalization group is then used to run between these scales 
to producing the resummed jet mass distribution. For simplicity, we extract these curves 
at fixed pt and y for the photon, although one could easily repeat this exercise integrated 
over some window of photon kinematics. 

6. Results 

To test our predictions, we compare our results to the output from PYTHIA 8 [12] for the 
LHC at £ C M = 8 TeV. We use the MSTW 2008 NLO PDFs [SI] both for event generation 
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Figure 8: The left panel shows a comparison of the leading order jet mass distirbution calculated 
with MCFM and the prediction from expanding the resummed result to leading order (SCET LO). 
This demonstrates that SCET reproduces the singularities of QCD at leading order. The right panel 
compares MCFM and SCET to pythia, with obvious disagreement, indicating the importance of 
resummation. 



and theory calculation. The event generation is made with the hadronization turned off, 
except for one curve on the right hand side of Figure where you can see that it affects the 
peak region. For the rest of plots, we keep the initial state and final state radiation, but 
we turn off hadronization, underlying event and multiple interactions. These additional 
effects are important, but we postpone their consideration to future work. The specific 
observable we use is outlined in Section 12.21 an d we test several values of R. The jet 
mass distribution is calculated and compared to pythia in a small window of the photon 
transverse momentum and rapidity (Apt, Ay) = (10 GeV, 0.1), centered around px = 500 
or 2000 GeV and y = 0, for various sizes of R. 

Since SCET contains all the physics in the small ttir regime, it should reproduce the 
singular structure of full QCD at leading order. We check this by comparing the fixed order 
expansion of SCET to the exact leading order distribution, using MCFM |62| . in the far 
singular region. The result is shown in Figure [U which shows very good agreement down 
to very small m,R. We also compare MCFM with pythia, and there is no region of uir in 
which the QCD NLO calculation agrees, which suggests the necessity of resummation. 

We calculate the jet mass distribution at full next-todeading log (NLL) precision, which 
involves the twodoop cusp anomalous dimension and the one loop anomalous dimensions of 
the hard, jet and soft functions. In fact, all ingredients are known for NNLL resummation, 
except for the contribution of non- global logarithms. We estimate the effect of these missing 
NGLs in Section [71 We denote the partial next-to-next-todeading log resummed result as 
NNLLp, which has the three loop cusp anomalous dimension and the twodoop anomalous 
dimensions of the hard, jet and full soft functions. We have only included those parts of 
the 0(a 2 s ) soft anomalous dimension that can be extracted from the RG invariance of the 
cross section, which does not include the angle-dependent pieces. 

Figure [9] shows the jet mass distributions for the px = 500 GeV events at different 
orders of precision. For this pr, we compute the jet mass distribution with a maximum 
cone size of R = 0.5 to avoid the jet being contaminated by the beam. We normalize the 
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Figure 9: Comparison of the jet mass distributions with different orders of precision (curves) 
for px = 500 GeV and R = 0.5 to pythia (histograms). The right plot shows the effects of 
hadronization in PYTHIA as compared to rcsummcd distributions. 



SCET distributions with the NLO QCD cross section as determined by pythia 8, scaled by 
the QCD NLO .fT-factor for the cross-section with the jet rapidity restriction (K ~ 0.8 for a 
500 GeV photon), which we compute with our own code. We do not match to the exact LO 
QCD calculation which would account the power corrections in the tail. In the peak region, 
where most of the events lie, these power corrections are small. The scale uncertainties for 
the resummed result include variation of the factorization scale fj,f = fx, the hard scale //^, 
the jet scale fj,j, and the soft scales /i in and // ut> Figure [Till shows the uncertainty bands for 
separate variation of the scales between ^/Xj < \i < 2//j for % = /, h,j, s- m , s out , for px = 500 
GeV and R = 0.5 jets. Additional comparisons for px = 2 TeV and R = 0.4 are shown 
in Figures [IT] and [12] The higher the transverse momentum of the photon, the closer to 
threshold, so we expect that threshold resummation will be more effective in this case. 

7. The Role of Non-Global Logarithms 

As mentioned in Section 2] although we were able to refactorize the soft function into a 
soft-collinear part, whose natural scale is associated with the soft modes within the jet, 
the remainder soft function still depended on multiple scales S r = S r (ki n , k out , //). Thus 
we cannot guarantee that all large logarithms in the jet mass distribution are resummed. 
The residual dependence of the remainder soft function on two scales is the problem of 
non-global structure. In the absence of a complete understanding of non-global structure, 
and how the non- global logs (NGLs) might be resummed, we will content ourselves with 
an estimate of how non-global structure might affect the jet mass distribution. We start 
by drawing on the lessons learned when considering e + e~ —> dijets, where there have been 
several studies [ISIIMIESIIISIET] . 

A two- loop calculation of the soft function for hemispherical jets was performed in [49 1 
[57] and for cone or anti-fej 1 jets, with an out-of-jet veto, in [25]. An intriguing observation 
from these studies is that in both examples, the NGLs arose from combinations of various 
logarithms of a single scale, coming from integrals in separate phase space regions. To 
give a specific example, in the calculation of the cumulative doubly differential dijet mass 
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Figure 10: Scale uncertainties for pr = 500 GeV and R = 0.5. The red bands are NLL and the 
blue bands are NNLL p . The PYTHIA result is the histogram shown in black. 



distribution for hemispherical jets, the details of which can be found in 
that the leading non-global logarithm is of the form 



it was shown 



3 A E R 



(7.1) 



This is the coefficient of (j^) in the soft function, with El and Er cumulant variables 
corresponding to integrals over Izl = n • Px L and Icr = n ■ Px R , respectively, where Px L 
and Px R are the total momentum of the soft radiation propagating in the left or right 
hemispheres and n = (l,n) denotes the thrust axis. This double logarithm came from the 
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sum of three contributions from the two gluon real emission graphs off the two Wilson lines 

CfCa In 2 — — both gluons in left hemisphere 

3 fi 

CfCa In 2 — — both gluon in right hemisphere • (7-2) 

3 ji 

H CfCa In — — In — — one gluon in each hemisphere 

3 [i [i 

One therefore expects a similar sum of phase space regions to produce the leading non- 
global logarithm for the soft function we need in this paper. 
Now, we have already refactorized the soft function as 

fx) = S(k in ,fi)S r (k in , k on t,fi) , (7.3) 

with the regional soft function S(ki n , fj,) defined to have all the radiation going into the jet. 
Thus we expect a double logarithm similar to the first term in Eq. (|7.2|) to contribute at 
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(7.4) 



where r NGL might be — ^-CfCa, but is yet unknown. Since this term has /U-dependence, 
it will contribute to the anomalous dimension of S. Moreover, it will contribute a term 
with In fi in the anomalous dimension. This is unusual since, for a global observable, one 
normally expects that all of the In [i dependence in the anomalous dimension is proportional 
to 

Tcusp- We therefore expect similar non-cusp In /i dependence in the anomalous dimension 
of the regional soft function for the jet mass distribution as well. 

With this educated guess, we expect that the two-loop regional soft function anomalous 
dimension for the annihilation channel has the form 

'1+/3 2 Kb 
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Utt/ 
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(7.5) 



and similarly for the Compton channel, with Cp 
two-loop cusp anomalous dimension: 
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and T NGL and 

^ngl are un k n0 wn. By RG invar iance, this implies that the anomalous 
dimension of the residual soft function must be 
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where 7 S is the two- loop direct-photon soft function anomalous dimension from [27] : 
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The terms proportional to T± and 7 S in Eq. (|7.5|) and Eq. (|7.7|) . and also a known term 
proportional to the 3-loop cusp anomalous dimension, account for the global logarithms at 
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Figure 13: Estimation of the effect of leading non-global log resummation. Note that it affects 
significantly only the peak region of the distribution. 

NNLL. The r NGL term parameterizing the unknown leading NGL. The 7 NGL anomalous 
dimension accounts for sub-leading NGLs, and the remaining non-global structure in the 
finite part of the two- loop soft function has effects which are formally N 3 LL. 

Without performing a two-loop calculation, the expressions for r NGL and 7 NGL are 
unknown. We estimate the effect of the missing terms by varying T NGL between ±100, 
since 100 ~ ACfT\ ~ ^-CfCa- The result of performing such a variation is shown in 
Figure [P51 The band in this plot shows a reasonable expectation of the improvement one 
could expect if the leading non-global logarithm could be resummed. We see that the NGL 
only affects significantly the distribution in the peak region. 

8. Conclusions 

In this paper we have calculated the distribution of jet mass at a hadron collider, for events 
in which the jet recoils against a hard photon. Our calculation includes resummation at 
the next-to-leading logarithmic level (NLL) which involves both final state radiation of the 
jet and initial state radiation of the colliding partons. We also have included resummation 
of all the global logarithmic terms at next-to-next-to-leading logarithmic level producing a 
distribution we label NNLL P ("p" for partial, since the non-global logarithmic terms have 
not been resummed). Our approach is based on expanding around the threshold limit, 
where the photon has very large momentum. By demanding the photon be hard, we force 
the hadronic final state to be that of a single jet for which a simple factorization formula 
exists. Our result is differential in the jet and photon rapidities and transverse momenta, 
although we assume the photon and jet have equal and opposite transverse momenta, which 
is true at leading power. 

We have compared our theoretical calculation to the output from pythia, and find 
very good agreement. Although pythia is only formally accurate to the leading-logarithmic 
level, it has elements of subleading logarithms, and is expected to be in good agreement 
with collider data. Since we are able to include color coherence effects, from initial state 
radiation, as well as subleading logarithms, our calculation is more precise than pythia, 
as far as perturbative QCD is concerned. Since PYTHIA also includes hadronization it is 
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more likely to describe the data well for very small jet mass. It would be interesting to 
compare both directly to collider data from the LHC when it becomes available. 

The two scales in our distribution, the photon px and the jet mass rriR lead to non- 
global logarithms. This is because in our approach there are two thresholds, itir — > and 
mx — > 0, where mx is the mass of everything-but-the-photon. In particular, there are two 
soft scales, h m associated with soft radiation in the jet, related to iur, and k out associated 
with radiation out of the jet, which can be related to m x — m 2 R . Non-global logarithms of 
the ratio of these scales make resummation difficult beyond the NLL level. An alternative 
calculation would be to impose a jet veto at a scale similar to the soft scale. Then the 
non-global structure would reduce to a single number, which one might estimate or argue 
to be small. However, for an inclusive jet mass calculation, it seems impossible to avoid 
non-global structure and therefore resumming non-global logarithms is necessary for NNLL 
resummation. 

Although we have not been able to resum the non-global logs in this paper, we found 
that if one ignores them completely, by choosing a single scale for the entire soft function, 
the distribution is completely wrong. This is because, using SCET alone, the two soft scales 
are not distinguished. Instead, we observe that contributions to the soft function coming 
from radiation going entirely into the jet can be consistently factorizes off. There are many 
ways to refactorize the soft function, and we make one particular choice, preserving the 
color structure associated with the hard directions. We give an operator definition to this 
regional soft function and calculate it at one loop. The regional soft function has a natural 
scale associated with soft/n j-collinear modes. Allowing ourselves to pick two different 
soft scales, which we do numerically, gives a result which is in very good agreement with 
pythia, in the region that our calculation can be trusted. 

There are many directions in which this work can be continued. It would be interesting 
to calculate the regional soft function at two loops, to see its non-global structure. One 
might hope that, since we expect an anomalous dimension for this soft function to have 
H dependence related to the leading non-global logarithm, that with further insight non- 
global logarithms could be resummed. Then one could produce an NNLL calculation of 
jet mass. Other related applications would include a calculation of the jet mass in dijet 
events, for which there is already data |63j and the one-loop hard function and anomalous 
dimensions have already been prepared [M], or a calculation of jet mass in direct photon 
events including a jet veto. With a jet veto one can force single-jet kinematics well away 
from threshold and the size of the non-global structure can be controlled. However, the 
calculation would be significantly more complicated than what we have done here. In 
addition, it would be interesting to pursue the consideration of the jet mass distribution 
using different jet algorithms. 
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